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The Letter by Krstulovic and Brachet (KB) [l[ ad- 
dresses an important problem in kinetics of Bose-Einstein 
condensation (BEC) in a weakly interacting Bose gas 
(WIBG) . By means of a numeric simulation of the trun- 
cated Gross-Pitaevskii equation (TGPE) (the high- wave- 
number harmonics are cut off at k ma , x ), the authors ob- 
serve a peculiar relaxation picture: starting with a su- 
perfluid vortex state, late-time evolution of the energy 
distribution towards thermal equilibrium takes on the 
form of a front at a characteristic wave number k c (t) (at 
which the energy is concentrated) with an abrupt trun- 
cation at the wave numbers fc max > k > k c [2[ . The front 
propagates toward higher k's at an ever-decreasing rate 
k c (t) leaving a quasithermalized distribution in its wake 
at k < k c . The Letter puts forward "a new mechanism 
of thermalization" suggesting that "a bottleneck delays 
the final thermalization when large dispersive effects are 
present at truncation wave number and produces an ef- 
fective self-truncation." We point out that the physical 
mechanism responsible for the numeric results is in fact 
quite different [3[: (i) in contrast to the bottleneck pro- 
posed in the delay in thermalization is due to the 
smooth increase of the typical kinetic time Tkm(^) with 
the wave number k, namely Tkin(fc) oc k, correspondingly 
(ii) the effective self-truncation is a consequence of energy 
conservation. The observed physics is well understood as 
a part of the general relaxation scenario of a strongly 
nonequilibrium WIBG developed by Svistunov [3|, de- 
scribing the underlying kinetic mechanisms, the form of 
the distribution, and its time dependence as k c (t) oc t x l A , 
which the numerics [Fig. 4(e) of Ref. [l|] agree with. 

The general relaxation scenario Q is nontrivial start- 
ing with an explosive wave towards k = 0, which popu- 
lates long-wavelength harmonics and leads to the forma- 
tion of the quasicondensate (QC) — the superfluid tur- 
bulence state, i.e. a tangle of quantized vortices — and a 
subsequent back wave propagating towards large fc's with 
a quasiequilibrium distribution formed behind the front. 
The simulation [l[ corresponds to the late-time stage of 
this scenario, in which most of the particles are already 
in the QC (represented in [l[ by a decaying Taylor-Green 


vortex) . This regime is described by the kinetic equation 
(KE) for the mode occupation numbers rik, 


rik = A r Coll ( [n k ),k\ 


(1) 


where Collo oc k rtf. is the collision integral and Nq is the 
number of QC particles, the specifics of QC dynamics be- 
ing irrelevant [4j . Nq is close to the total particle number 
and can be considered constant. Physically, the KE de- 
scribes the dominant three-wave collision processes with 
the fourth wave being the QC iVo- Scale invariance of 
the KE along with the conservation of the total energy 
E oc J k 4 dk rife dictates that the solution describing the 
thermalization at high fc's takes on a self-similar form 
(Eq. (4.4) in Ref. Q): 


n k = k c (t)- 5 f(k/k c (t))- 


(2) 


Here f(x) — 0, x ^> 1 enforcing truncation at k c to 
ensure convergence of E, and f(x) oc 1/x 2 , x <C 1 corre- 
sponds to the Gibbs distribution for k <C fc c - The evolu- 
tion of k c (t) follows from the KE (JTJ> yielding k c (t) oc t 1 / 4 . 

Eqs. (fT]), © prescribe that kinetics get gradually 
slower with k, 7ki n (&) oc k. That is the fundamental 
reason for the deceleration of the front propagation at 
k c (t) and the thermal equilibration of the modes in its 
wake. The solution $2$ describes a drift of energy [con- 
trasted with an energy cascade, in which T\,i n (k) vanishes 
towards k — > oo) [3|, whereas a bottleneck, i.e. an abrupt 
drop in kinetic efficiency at some k-scale (typically due to 
a change of the kinetic mechanism], does not take place. 

KB also raise an interesting question of whether the 
ultimate thermalization is completely inhibited in the 
limit of an arbitrarily high spatial resolution due to the 
ever-slowing-down nature of the process at high fc's. In 
real systems, the classical-field regime of the GPE breaks 
down at the wave numbers for which n k ~ 1 and the ki- 
netics are dominated by the quantum spontaneous (as op- 
posed to stimulated at S> 1) scattering processes 0]. 
The spontaneous scattering provides an efficient mecha- 
nism for the equilibration at high wave numbers. 
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TGPE due to irrelevance of the modes with k ~ k n 
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